Report for Model 1

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.01       1.02
## dh0           1.42       2.05
## dh1           1.13       1.36
## dl0           1.51       2.22
## dl1           1.06       1.16
## 
## Multivariate psrf
## 
## 1.61

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1343.40738 22794.8154
DIC3 1218.27290 22823.8719
PWAIC 63.37221 196.0088
WAIC 1258.22345 22835.7281

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.00
## LDevsum       1.00       1.00
## dh0           1.84       3.05
## dh1           1.86       2.98
## dl0           1.28       1.71
## dl1           1.09       1.26
## 
## Multivariate psrf
## 
## 1.71

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1353.59804 33770.3243
DIC3 1230.04462 33889.3234
PWAIC 84.87555 167.7255
WAIC 1295.26122 33894.1349